Tunneling mechanism of light transmission through metallic films 
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A mechanism of light transmission through metallic films is proposed, assisted by tunnelling be- 
tween resonating buried dielectric inclusions. This is illustrated by arrays of Si spheres embedded in 
Ag. Strong transmission peaks are observed near the Mie resonances of the spheres. The interaction 
among various planes of spheres and interference effects between these resonances and the surface 
plasmons of Ag lead to mixing and splitting of the resonances. Transmission is proved to be limited 
only by absorption. For small spheres, the effective dielectric constant can be tuned to values close 
to unity and a method is proposed to turn the resulting materials invisible. 

PACS numbers: 78.66.Bz,42.25.Fx,73.20.Mf,42.79.Ci 



Substantial efforts have been placed on studying light 
transmission through perforated metallic films since the 
discovery by Ebbesen et at fl| of extraordinary transmis- 
sion through sub-wavelength periodic hole arrays, lead- 
ing to such remarkable effects as giant funnelling of mi- 
crowave radiation through metallic micrometer slits ,2]. 
Two distinct mechanisms have been identified that con- 
tribute to enhance light transmission: (i) dynamical scat- 
tering assisted by surface plasmons (SPs) and (ii) reso- 
nant transmission by coupling to propagating modes, (i) 
Isolated sub-wavelength holes can only support evanes- 
cent modes, but dynamical diffraction assisted by SPs in 
periodic hole arrays lead to transmission maxima that 
can beat the severe 1/ dependence on wavelength A of 
the transmission derived by Bethe for single apertures 
0, although there is still some controversy regarding 
the relative relevance of SPs and dynamical diffraction 
BHHIl- (ii) Slits differ from holes in that the former, 
however small, support at least one transmission mode 
that can couple resonantly to external light to transmit 
several orders of magnitude more light than impinging 
directly on the slit aperture 2] . 

In this Letter we investigate yet another mechanism of 
enhanced transmission based upon hopping between di- 
electric inclusions that can sustain localized resonances 
inside metallic films. We illustrate this concept by 
studying light transmission through metallic Ag films 
that contain spherical Si inclusions arranged in square- 
lattice layers, which belongs to the family of periodic 
arrays of nanoresonators coupled via spatially homoge- 
neous waveguides Q . The transmission is shown to take 
large values near the Mie resonances of the spheres, and 
it is only limited by metal absorption. 

We first consider a single layer of inclusions [Fig.[T]Ja)]. 
Maxwell's equations are solved rigourously for this geom- 



etry using a layer-by-layer version of the KKR multiple- 
scattering method [9], with lattice sums performed in 
real space within the metal. The complex, frequency- 
dependent dielectric functions of Ag and Si are taken 
from optical data [lol |. 

The resonant modes of the Si spheres can be visual- 
ized through their Mie scattering coefficients inside bulk 
metal [IL] , showing prominent features in the visible and 
near-IR domains for a sphere radius R = 75 nm [Fig. 
[IJb)]. The corresponding transmittance spectrum ex- 
hibits peaks near those resonances [Fig. [H^c) , at normal 
incidence] that can exceed a transmission T — 50% at 
wavelengths near A = 975 nm, well above the T = 13% 
value expected for a homogeneous Ag film of thickness 
2b = 15 nm (i.e., the combined thickness of the thin bury- 
ing layers at either sides of the spheres). Smaller separa- 
tions between the sphere and the film boundaries lead to 
even higher transmission (e.g., T = 70% for b — 2.5 nm). 
The transmission peaks are red-shifted with respect to 
the resonances in bulk Ag as a result of interaction be- 
tween sphere modes and vacuum modes. However, the 
relevant multipoles of the inclusions and the polarization 
symmetry associated to those peaks can be well ascribed 
by comparing Figs.[ljb) and (c). 

The interaction of spherical inclusions along the plane 
can be important at small separations between the 
spheres. This is illustrated in Fig. [2l which shows a con- 
tour plot of the transmittance T as a function of lattice 
constant a and wavelength. T has been multiplied by 

/ttE? to obtain transmission cross sections per sphere 
The El and E2 Mie resonances of Fig. [IJb) show up 
as two lines of transmission maxima relatively indepen- 
dent of a. However, these two modes interact and repel 
each other when the spheres are near touching each other 
(bottom part of the figure). 
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FIG. 1: (a) Transversal section of a Ag film containing a 
buried square planar array of 150-nm Si spheres in square 
lattice configuration, (b) Square modulus of the Mie scatter- 
ing coefficients of single 150-nm Si spheres inside bulk Ag for 
different multipoles {I — 1,2,3) and polarization modes [elec- 
tric (E) and magnetic (M)] as a function of light wavelength 
in vacuum, (c) Absolute transmittance through the film of 
(a) under normal incidence for a lattice constant a — 160 nm. 
The insets show maps of the electric field strength in the slab 
in logarithmic scale for circularly polarized incident light at 
the peak wavelengths. 
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FIG. 2: Transmittance through the film described in Fig.ma) 
under normal incidence as a function of wavelength and lattice 
constant. The transmittance is normalized to the projected 
area of the spheres. The solid curves are given by the condi- 
tion that the momentum of the planar surface plasmon of Ag 
matches some reciprocal lattice vector of the sphere lattice. 
The right inset is a blow up of an interference region and the 
left inset shows the cross section along the segment AB. 



The SPs of the planar Ag interfaces cannot couple 
directly to external propagating modes. However, the 
presence of Si inclusions provides a source of momentum 
that makes this coupling possible. The coupling increases 
when the interaction between lattice sites (spheres) via 
SPs is constructive, which at normal incidence (all lat- 
tice sites are in phase) leads to the kinematical condi- 
tion that SP momenta that match reciprocal lattice vec- 
tors of the square lattice of spheres are favored, that is, 
AV + w? = a\J CAg/ (cAg + 1) [H, where the surface 
plasmon dispersion relation has been invoked. The solid 
curves in Fig. [2] are obtained from this equation for differ- 
ent combinations of integers n and m and they would lie 
near regions of transmission enhancement in hole arrays 
[H , but that transmission mechanism is marginal for our 
buried structures. Instead, a rich transmission structure 
shows up when the noted kinematical condition (solid 
curves) is at the same wavelength as a Mie resonance [see 
right inset of Fig. [2]: the discrete, localized Mie modes 
couple to the continuum of SPs to yield a transmission 
spectrum exhibiting a point of zero transmission, as it is 
well known from Fano's theory (see left inset of Fig. 
[2]). 

Several layers of spheres can interact in a single film 
producing splitting and shifting of the Mie resonances, 
as shown in Fig. [31 The main 975-nm feature of Fig. 
[IJc) is subdivided here into a series of N peaks in a film 
containing N layers. The coupling between Mie modes 
at consecutive layers is the same as in previously intro- 
duced photonic molecules that arise from the interaction 
of Mie modes in neighboring dielectric particles [l^ . The 
tight-binding approach used in that context becomes a 
natural tool to understand the hopping of light between 
contiguous spheres in our case. In this spirit, an anal- 
ogy with electron transmission between localized states 
in quantum systems can be established by considering 
a simple tight-binding model [14] in which a resonant 
mode trapped in each layer is associated to a state |i) 
(with i ~ 1,2, ...,N). The dynamics is governed by the 
Hamiltonian 

H = {u;b + i^sm{l\ + \N){N\)+ {LJb + iJ:b)\i){i\ 

i=2 

i=N-l 

+ E {V\i){i + l\ + h.c.), 



where Wf, and Sf, are the energy and width of the reso- 
nances, and V is the hopping amplitude, which expresses 
the magnitude of the interaction between spheres. For 
internal layers the width Ef, accounts for dissipation into 
bulk metal. The states |1) and |A^) interact with the con- 
tinuum of light modes on the left and right hand sides 
of the film, respectively, described by frequency depen- 
dent self-energy Es(tj). The transmittance in the model 
is defined by the squared absolute value of the transition 



3 



A' layers 



b=1.5 nm-i '^d=S nm 




0.0 
850 



1 1-^ 


A 















O 
03 




850 



900 950 

Wavelength (nm) 



1000 



900 950 

Wavelength (nm) 



FIG. 4: Absolute transmittance under the same conditions 
as in Fig. [3l but setting the imaginary part of the dielectric 
constant of Ag to zero. 



FIG. 3: Transmittance through Ag films containing different 
numbers of buried planar arrays of Si spheres of radius R — 
75 nm under normal incidence conditions. The planar unit cell 
is a square of side a = 250 nm, which prevents interaction 
between spheres within each plane. Solid curves are exact 
solutions of Maxwell's equations. Broken curves are obtained 
from the tight-binding model explained in the text [Eq. ([l])]. 

amplitude and can be written as 

T-4|S,H.gi.Ar(Lc> + iO+)|2, (1) 

where (7i_7v(^) = {mz~H)~^\N). The parameters of the 
Hamiltonian have been chosen to produce the best fit to 
the transmittance obtained from numerical solution of 
Maxwell's equations: S;, = — 0.002a;h, V = O.OScjf,, and 
a linear function coupling to vacuum given by 'Ss{'-^) = 
—0. 005^6 — 0.03(Li; — ujb)- Note that radiative couphng 
via non- vanishing Im{y} can produce changes in the line 
shapes [l5j] that lie beyond the scope of this work. 

The results derived from Eq. ([1]) are shown in Fig. [3] 
as dashed curves. Characteristic features of the exact 
transmittance (solid curves) are well reproduced, includ- 
ing peak splitting, peak asymmetry coming from the en- 
ergy dependence of Ss , and narrowing of the resonances 
with increasing number of layers. This model corrobo- 
rates that the main mechanism producing transmission 
is tunnelling between electromagnetic modes trapped in 
each layer. 

The dielectric function of Ag near the main resonance 
of Fig.[3]is essentially negative (eAg = — 49-f 0.6i), so that 
coupling among the spheres and between the spheres and 
the film planar surfaces is mediated by evanescent modes. 
However, the small imaginary part of eAg plays a leading 
role in reducing the absolute transmittance (Im{esi} <C 1 



fld^). This is made clear in Fig. [H which has been calcu- 
lated under the same conditions as in Fig. [3] but setting 
the imaginary part of eAg to a positive infinitesimal. All 
resonances in the figure reach a maximum of 100% trans- 
mission, a fact that can be rigorously explained following 
simple arguments based upon the change by tt of the 
phase of reflected and transmitted resonant components 
of the field across the resonance [16]. 

A Mie cavity mode inside a lossless bulk metal must 
have infinite lifetime because light can neither scape 
through nor be absorbed by the metal. Therefore, the 
width of the resonance in the lossless film for iV = 1 
(6.72 nm FWHM, as obtained from Fig. |4]) has to origi- 
nate in the coupling of the El mode with the continuum 
of free states on either side of the film (leakage into the 
vacuum). For multilayer systems (e.g., N — 2), as a 
consequence of trace preservation in interacting systems 
as compared to the non-interacting ones, this broadening 
has to remain the same (indeed, the sum of the FWHM in 
the two peaks for = 2 is 6.72 nm as well). This leakage 
contribution to peak broadening adds up to the intrinsic 
width of the Mie resonance in the lossy bulk metal [4.53 
nm FWHM as obtained from the El peak in Fig. [Ijb)], 
yielding a total width of 11.25 nm, which compares well 
with the value of 11.83 nm observed for the actual lossy 
metallic film {N = 1 solid curve in Fig. [3]). 

Therefore, our metallic film would be invisible to nor- 
mally incident radiation at specific resonant wavelengths 
if the metal was free of absorption. The question arises 
whether the same mechanism of hopping between dielec- 
tric inclusions in metals can be utilized to fabricate solid 
structures that are invisible to light incident with any 
angle at some wavelength, or in other words, whether 
a medium exists exhibiting effective dielectric function 



4 



■3 

> 





Re{frff} 




Im{£'|,ff} for the 
pasive material 


\ 

Im{e,ff} for 
active mate 


the 
■ial 

Ag — 


m 


Si particles^ v,^^ 
doped with 
active atoms 



720 



740 760 

Wavelength (nm) 



780 



FIG. 5: Effective dielectric function of an fee array of Si spher- 
ical inclusions in Ag with a filling fraction of Si of 50% (solid 
curves). The effective dielectric function is modified when 
Si is doped with active centers (broken curves), described by 
adding a negative imaginary part (—0.2521) to the its dielec- 
tric function. The resulting material is invisible (eeff ~ 1) at 
a wavelength of 735 nm. 



eoff = 1 (and magnetic response ficS = 1 as well, of 
course). To explore this possibility, we have represented 
Eoff in Fig.[5]for a composite material containing 50% of Si 
inside spherical inclusions, which are arranged within Ag 
in an fee lattice and separated by distances much shorter 
than the wavel eng th, eos is obtained using Maxwell- 
Garnett theory [17| , and we have found qualitatively sim- 
ilar results from rigorous solution of Maxwell's equations 
for the reflectivity of a surface of the composite material 
as compared to Fresnel equations for an equivalent ho- 
mogeneous medium. There is a wavelength at which the 
real part of ees is 1, but its imaginary part takes a small 
but non-negligible value. 

We propose to dope Si in such a structure with 
optically-active atoms or molecules that are capable of 
sustaining an inversion of population under excitation 
by light of a different wavelength or by electric discharge. 
The resulting lasing activity can be well represented by 
a negative imaginary part contributing to the dielectric 
function of doped Si [18|, whose value can be controlled 
by the doping dose or, in 4-level active atoms, also via 
the pumping intensity. In particular, when this imagi- 
nary contribution to egi is taken as — 0.252i then ecff = 1 
at 735 nm (the imaginary part of eos cancels out exactly), 
and the material becomes invisible at that wavelength. 



In conclusion, we have shown that light transmission 
through metallic films can be enhanced by coupling to 
resonances localized in dielectric inclusions, paving the 
way towards driving light deep inside metals. Light 
circuits based upon this mechanism can be envisioned. 
When the inclusions are small compared to the wave- 
length, the resulting metamaterial can be made invisible 
(i.e., eeff = 1) if metal absorption is compensated by dop- 
ing the dielectric with active atoms under inversion pop- 
ulation conditions. We have proposed a specific design 
for such invisible material at wavelengths in the visible, 
but structures insensitive to microwaves should be also 
realizable due to smaller absorption in this range and 
using reradiating active elements in the millimeter scale. 



T. W. Ebbesen et al, Nature 391, 667 (1998); H. F. 
Ghaemi et al, Phys. Rev. B 58, 6779 (1998). 
Y. Takakura, Phys. Rev. Lett. 86, 5601 (2001); F. Yang 
and J. R. Sambles, Phys. Rev. Lett. 89, 63901 (2002). 
H. A. Bethe, Phys. Rev. 66, 163 (1944). 
M. M. J. Treacy, Appl. Phys. Lett. 75, 606 (1999); Phys. 
Rev. B 66, 195105 (2002). 

E. Popov et al, Phys. Rev. B 62, 16100 (2000). 
L. Martm-Moreno et al, Phys. Rev. Lett. 86, 1114 
(2001). 

M. Sarrazin et al, Phys. Rev. B 67, 085415 (2003). 
S. G. Tikhodeev et al, Phys. Rev. B 66, 045102 (2002); 
A. Christ et al, Phys. Rev. Lett. 91, 183901 (2003); T. 
Zentgraf et al, Phys. Rev. Lett. 93, 243901 (2004). 
N. Stefanou et al., Comput. Phys. Commun. 113, 49 
(1998); 132, 189 (2000). 

E. D. Palik, Handbook of Optical Constants of Solids 
(Academic Press, New York, 1985). 

J. D. Jackson, Classical Electrodynamics (Wiley, New 
York, 1975). 

U. Fano, Phys. Rev. 124, 1866 (1961). 
T. Mukaiyama et al., Phys. Rev. Lett. 82, 4623 (1999). 
E. N. Economou, Creen's Functions m Quantum Physics 
(Springer- Verlag, Berlin, 1983). 

C. Ropers et al, Phys. Rev. Lett. 94, 113901 (2005). 
A. G. Borisov et al, Phys. Rev. B 71, 075408 (2005). 
J. C. Maxwell-Garnett, Philos. Trans. R. Soc. London A 
203, 385 (1904); 205, 237 (1906). 
[18] See V. I. Kopp et al, Phys. Rev. Lett. 86, 1753 (2001). 
We overlook the change in Re{esi} here, as it only pro- 
duces a shift of the wavelength of invisibility. 



[1 

[2: 

[3] 
[4] 

[6: 

[7 

[9: 

[10 

[11 

[12 

[13 
[14 

[15 
[16 
[17 



